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Abstract. The structure of the high galactic latitude Cometary Globule 12 (CG 12) has been investigated by means of radio 
molecular line observations. Detailed, high signal to noise ratio maps in C ls O (1-0), C 18 (2-1) and molecules tracing high 
density gas, CS (3-2), DCO+ (2-1) and H 13 CO + (1-0), are presented. The C ls O line emission is distributed in a 10' long 
North-South elongated lane with two strong maxima, CG 12-N(orth) and CG 12-S(outh). In CG 12-S the high density tracers 
delineate a compact core, DCO + core, which is offset by 15" from the C 18 maximum. The observed strong C 18 emission 
traces the surface of the DCO + core or a separate, adjacent cloud component. The driving source of the collimated molecular 
outflow detected by White (1993) is located in the DCO + core. The C 18 lines in CG 12-S have low intensity wings possibly 
caused by the outflow.The emission in high density tracers is weak in CG 12-N and especially the H 13 CO + , DCO + and 
N2H + lines are +0.5 kms -1 offset in velocity with respect to the C ls O lines. Evidence is presented that the molecular gas 
is highly depleted. The observed strong C 18 emission towards CG 12-N originates in the envelope of this depleted cloud 
component or in a separate entity seen in the same line of sight. The C ls O lines in CG 12 were analyzed using Positive Matrix 
Factorization, PMF The shape and the spatial distribution of the individual PMF factors fitted separately to the C ls O (1-0) 
and (2-1) transitions were consistent with each other. The results indicate a complex velocity and line excitation structure 
in the cloud. Besides separate cloud velocity components the C ls O line shapes and intensities are influenced by excitation 
temperature variations caused by e.g, the molecular outflow or by molecular depletion. Assuming a distance of 630 pc the size 
of the CG 12 compact head, 1.1 pc by 1.8 pc, and the C 18 mass larger than 100 Mq are comparable to those of other nearby 
low/intermediate mass star formation regions. 

Key words, clouds - ISM molecules - ISM: structure - radio lines - ISM: individual objects: CG 12, NGC5367 



1. Introduction 

Herschel (1847) noted that a 10 th mag. star, now known as 
h4636 or CoD -39° 8581, is a binary with a separation of 
3'/7. In the optical, h4636 illuminates the bright reflection neb- 
ula NGC 5367. |Hawarden & Brand (1976)| showed that it lies 
in the head of an impressive Cometary Globule 12, CG 12, 
with a tail stretching about one degree to the SE. With a 
galactic latitude of 21° and at the distance of 630 pc esti- 
mated by |Williams et al. (1977)1 CG 12 lies more than 200 pc 
above the plane. CG 12 has an associated low/intermediate 
mass stellar cluster which has at least 9 members (Williams 
et al. 1977). The clouds cometary structure could be due to 
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the passage of a supernova blast wave. Curiously, the cometary 
tail stretches towards the Galactic plane which would place 
the putative supernova even farther away from the Galactic 
plane than the globule. According to |Maheswar et al (1 996) 
the head of CG 12 is pointing towards the centre of an HI 
shell. Such a shell is, however, not readily evident in the 
whole sky HI survey ( Kal berla et al. 2 005 ) which merges the 
northern Leiden/Dwingeloo Survey ( Hartman & Burt on 1997 \ 
and the southern Instituto Argentino de Radioastronomia 
Survey dArnal et al. 2000> . The CG 12 cloud cometary shape 
is also seen in the IRAS surface emission (Oden wald 19881 . 

White (1993) mapped the region around and south of h4636 
in 12 CO (2-1) and C ls O (2-1) lines with a spatial resolution 
of 22". He found a small size C ls O core near the binary, and 
a molecular outflow with a centre close to the binary system. 
Further large scale 12 CO , 13 CO and C 18 observations with 
2.'7 resolution are presented in Yone kura et al. (1999)| CG 12 
contains two compact 1.2mm continuum sources, one in the 



2 



L. K. Haikala and M. Olberg: The structure of the cometary globule CG 12 




13 h 57 m 40" 30" 20" -3 -4 

R.A. (J2Q00) R.A. offset (arcmin) 



Fig. 1. Maps of integrated emission of C ls O (1-0) and C ls O (2-1) / T* A dv from -8 kms" 1 to -4 kms" 1 in CG12. The lowest 
contour and the contour increment are 0.3 Kkms -1 . The SEST half power beam sizes are shown in the upper corners of the 
panels. Offsets from the map centre position, 13 h 57 m 43!l, -39°58'43'/3 (J2000), are shown on the axes of the right panel. The 
scale tick in the lower right corner assumes a distance of 630 pc to the cloud. 



direction of the C 18 core detected by White (1993) (Reipurth 
et al. 1996) and another one two arcminutes north of it (Haikala 
2006, in preparation). The centres of the continuum sources 
were observed in C ls O (3-2) by Haik ala et al. (2006)| 

Near infrared (J, H and K) images and photometry of 
stars in CG12 are available in the 2MASS survey and 
|Santos et al. (1998)1 A deeper J, H and Ks imaging study is 
presented in Haikala (2006, in preparation). Far infrared emis- 
sion in CG 12 is dominated by a strong point source, IRAS 
13547-3944, near the binary h4636. 

CG 12/NGC 5367 is an intriguing object. It has the ap- 
pearance of a cometary globule, such as usually found in the 
outskirts of HII regions. The linear size of cometary glob- 
ules, like those in the Gum nebula, are however much smaller 
(e.g., Reipurth 1983). The linear extent of CG12, 10 pc, is four 
times larger than, e.g., the archetype object CG 1. A cluster of 
low/intermediate mass stellar cluster is associated with CG 12. 
Perhaps the most curious feature of CG 12 is its location over 
200 pc above the Galactic plane with no sign of other nearby 
dark clouds or star formation. CG 12 has not attracted much 
interest since the identification of the stellar cluster (Williams 



et al.1977) and the Hawarde n & Brand (1976)| cometary glob- 
ule paper. The subsequent papers have either concentrated on 
the Herbig AeBe binary, h4636 or the cloud has been included 
in various surveys. The 12 CO (2-1) and C 18 (2-1) observa- 
tions by White (1993), though detailed, cover only the very 
centre of the cloud. What has been missing is a detailed but at 
the same time extended study of the dense CG 12 molecular 
cloud in molecular transitions sensitive to the large scale struc- 
ture (cloud envelope) and to the detailed structure of the high 
density material (cores). 

In this paper we report the mapping of the head of 
CG12 in the C ls O (1-0) and (2-1) and in 13 CO (1-0) 
lines. The C 18 emission maxima were further mapped in 
DCO+(2-l), H 13 CO+(1-0) and CS (2-1) and (3-2). 
Further pointed observations with long integration time in 
CO (1-0) and (2-1) (and isotopologues), CS (2-1) and (3-2), 
C 34 S (2-1), H 13 CO+ (1-0), DCO+(2 - 1) and N 2 H+ (1-0) 
were made towards selected positions in the cloud. 

Observations, data reduction and calibration procedures are 
described in Sect. |2] and the observational results in Sect [3] 
The new results are compared with the optical and NIR im- 
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ages in Sect. 0] In Sect. [5] and Appendix |X] Positive Matrix 
Factorization is used to analyze the C 18 small scale struc- 
ture in the globule. In the discussion part in Sect. [6] the 
C 18 column densities and the cloud mass are derived and the 
observational and calculated results are summarized. The con- 
clusions are drawn in Sect.0 

2. Observations 

The observations were made during various observing runs 
with the Swedish-ESO-Submillimetre-Telescope, SEST, at the 
La Silla observatory, Chile. The SEST 3mm dual polariza- 
tion, single sideband (SSB), Schotky receiver was used for the 
13 CO (1-0) and CS (2-1) mapping observations. Rest of the ob- 
servations were conducted with the SEST 3 and 2 mm (SESIS) 
and 3 and 1 mm (IRAM) dual SiS SSB receivers. The SEST 
high resolution 2000 channel acousto-optical spectrometer 
(bandwidth 86 MHz, channel width 43 kHz) was split into two 
halves to measure two receivers simultaneously. At the ob- 
served wavelengths, 3mm, 2mm and 1mm, the 43 kHz chan- 
nel width corresponds to ~0.12 kms -1 , ~0.08 kins" 1 and 
^0.06 kms -1 , respectively. 

Frequency switching observing mode was used and a sec- 
ond order baseline was subtracted from the spectra after fold- 
ing. Calibration was achieved by the chopper wheel method. 
All the line temperatures in this paper, if not specially noted, 
are in the units of T\ , i.e. corrected to outside of the atmo- 
sphere but not for beam coupling. Typical values for the effec- 
tive SSB system temperatures outside the atmosphere ranged 
from 200 K to 350 K. Pointing was checked regularly in con- 
tinuum mode towards the nearby Centaurus A galaxy. Pointing 
accuracy is estimated to be better than 5". 

The observed molecular transitions, their frequencies, 
SEST half power beam width, HPBW, and the telescope main 
beam efficiency, r? m b, at these frequencies are given in Tabled 
For the CO observations only the C 18 is listed. At a distance 
of 630 pc the SEST HPBW at 219 GHz corresponds to 0.07 pc. 

The cloud was mapped simultaneously in C ls O (1-0) and 
C ls O (2-1) transitions with a spacing of 20" (514 positions) . 
The map centre was 13 h 57 m 43!l, -39°58'43"3 (J2000) which 
is —5" and ~ —10" in right ascension from the positions of the 
IRAS 13547-3944 point source and the binary h3626, respec- 
tively. The average rms of the spectra were 0.07 K and 0.09 K 
for C 18 (1-0) and (2-1), respectively. An approximately 8' 
by 20' area was mapped in 13 CO using 40" spacing (324 po- 
sitions). The C ls O maxima were mapped in CS (2-1), (3-2), 
DCO+(2-l) and H 13 CO+(1-0). Further long integration time 
pointed observations in CO, CS, H 13 CO+, DCO+ C 34 S and 
N2H + were made. 

3. Results 

3. 1. CO 

The observed distributions of C ls O (1-0) and C ls O (2-1) line 
emission towards CG12 are presented in Fig. Q(the JT^ dv 
in the velocity range —8 kms -1 to —4 kms" 1 ) and Fig. II II 
(C ls O channel map). In both figures the grey/colour scales and 



Table 1. Observed lines and telescope parameters 



Line 


v [GHz] 


HPBW ["] 




H la CO+ (1-0) 


86.754 


57 


0.75 


N 2 H+ (1-0) 


93.176 


54 




C 34 S (2-1) 


96.412 


54 




CS (2-1) 


97.271 


54 




C 18 O(l-0) 


109.782 


47 


0.70 


DCO+ (2-1) 


144.077 


34 




CS (3-2) 


145.904 


34 


0.66 


C 18 (2-1) 


219.560 


24 


0.50 



the contour levels are the same for the two C 18 transitions. 
The offsets from the map centre position in arc minutes are 
shown on the axes of the right panel in Fig.^ 

The bulk of the molecular material traced by C ls O emis- 
sion is distributed in a narrow North-South oriented lane with 
two prominent maxima. A less intense maximum is observed to 
the SW. Henceforth the three intensity maxima will be referred 
to as CG 12-N, CG 12-S and CG 12-SW. CG 12-S corresponds 
to the C 18 (2-1) maximum reported in |White (1993)1 

The morphology of the cloud is similar in the two 
C 18 transitions. However, there is a notable difference be- 
tween the C 18 maxima. The observed C ls O (1-0) emission 
is stronger than the C ls O (2-1) emission in CG 12-N and 
CG 12-SW whereas in CG 12-S the opposite is the case. In 
CG 12-N the C 18 (1-0) line remains stronger than the (2-1) 
line even when expressed in the main beam brightness temper- 
ature scale. 

Further details can be seen if Fig. which reveals that 
the structure of the cloud is not as simple as Fig. ^ suggests. 
The most noticeable features are the following: CG 12-N is 
elongated in the North South direction at velocities < —6.0 
kms -1 and in the East West direction at velocities > —5.4 
kms -1 . The position of the maximum C 18 (2-1) emission in 
CG 12-S moves from below the map centre position (marked 
with a cross in the Figure) at —7.0 kms -1 to a position 40" 
west of this position at velocity —5.9 kms -1 . An arc like fea- 
ture connects CG 12-S and CG 12-SW at velocities from —6.6 
kms -1 to —6.2 kms -1 . The arc is less pronounced in the 
C 18 (2-1) transition. 

The observed distribution of the 13 CO (1-0) line emis- 
sion towards CG12 is presented in Fig. |2] The extent of the 
C 18 mapping and the outline of the C 18 (2-1) emission is 
indicated in the overlay. The notable difference between the 
observed C 18 and 13 CO emission is that there is only one 
13 CO maximum which is offset to NE from CG 12-S. In par- 
ticular there is no indication of CG 12-N in the 13 CO map. 

12 CO, 13 CO and C ls O (1-0) and (2-1) (pointed, long inte- 
gration time) spectra in five selected positions in the cloud are 
shown in Fig. [3] The 12 CO spectra observed at all positions in 
the figure are strongly self-absorbed. Therefore the 12 CO line 
peak intensity, line half width and the line integral are not phys- 
ically meaningful. Line wing emission due to a molecular out- 
flow (White 1993) is seen in all 12 CO spectra, also in the direc- 
tion of CG 12-N which was not covered by the |White (1993)| 
observations. 
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Fig. 3. 12 CO , 13 CO and C 18 spectra observed towards five selected positions in CG 12. Offset from the map centre position 
is shown in arc seconds in the upper left corners of each panel. Line has been used for the (1-0) and histogram for the (2-1) 
transitions. The antenna temperature scale at left is for 12 CO and 13 CO spectra and the one on right for C ls O. The velocity tick 
indicates velocity of —6.5 kms -1 . The 12 CO and 13 CO spectra in position (—40", 160") are offset from the zero level for clarity 



3.2. High density tracers, mapping 

In the optically thin case the critical density of the C 18 (1— 

0) transition is ^650 cm -3 and approximately ten times 
this value for the (2-1) transition dRohlfs & Wilson 19991 . 
These C ls O lines are therefore excited already at low den- 
sities and their emission traces column density rather than 
number density. Therefore the C 18 intensity maxima were 
also mapped in the high density tracers (critical densi- 
ties larger than 10 5 cm- 3 ) CS (2-1), H 13 CO+(l - 0) and 
DCO + (2 — 1) using 20" spacing. A 4 by 4 point map with the 
same spacing was obtained in CG 12-S in the CS (3-2) line. 

Contour maps of the CS (2-1), CS (3-2), H 13 CO+ and 
DCO+ line integral in CG 12 superposed on the C 18 (2- 

1) emission (grey scale) are shown in Fig. |4] CG 12-N was 
not mapped in CS (3-2). The CS (2-1) emission peaks at 
the position of the C ls O (2-1) the maximum in CG 12-S. 
In other molecules the maximum is shifted to SE from 
the C 18 maximum. As DCO + traces high density gas the 
CG 12-S DCO+ maximum will be referred to as the DCO+ 
core. 



3.3. C 18 and high density tracers: pointed 
observations 

Pointed, long integration time CS (2-1), (3-2), C 34 S(2-1), 
H 13 CO+ (1-0), DCO+ (2-1) and N 2 H+ (1-0) spectra in the 
same positions as in Fig [5] are shown in Fig. [5] The integra- 
tion time of the CS (3-2) line at position (-40", -20") is only 
2.5 minutes and is therefore noisier than the two other CS (3-2) 
lines. N2H + (1-0) line consists of seven hyperfine components 
and only the unblended N 2 H+ (1-0) (Fi, F= 0, 1 -> 1, 2) com- 
ponent is shown in Fig. [5] 

The C 34 S line was observed in three positions. These spec- 
tra together with the CS (2-1) and (3-2) spectra observed at the 
same positions are shown in Fig. [6] 

The N2H+ spectra showing all the hyperfine components 
in the three observed positions are shown in Fig. Also the 
hyperfine component fits (assuming optically thin emission and 
one velocity component) and their residuals are shown. 
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Fig. 2. Map of integrated emission of 13 C0 (1-0): JT^ dvfrom 
—8 kms -1 to —4 kms -1 . The lowest contour and the contour 
increment is 1.2 K kins" 1 . The SEST half power beam size at 
the line frequency is shown in the upper corner. The extent of 
the C 18 mapping and the outline of the C 18 (2-1) emission 
is indicated in the overlay. 



3.3.1. CG 12-S 

Near the centre of the DCO+ core (position -20",— 20") 
the H 13 CO + , DCO + and N 2 H + lines are nearly symmet- 
ric and centered at velocity —6.5 kms -1 . The C 18 lines 
are skewed and the maximum intensity is redshifted from 
the DCO + (2 — 1) peak velocity. Going to the West (position 
-40", -20") the redshifted side of the C 18 lines becomes 
more intense. Low intensity wing-like emission is observed in 
the C 18 and CS (2-1) lines in the direction of the DCO+ core 
and position (0",0"). 

The N 2 H+ line in the position (-20",-20") can be fit with 
a single velocity component centred at —6.56 kms -1 (Fig.Q 
lower panel). Due to the 0.73 km s _1 line half width the hyper- 
fine components are blended and the hyperfine structure is only 
marginally resolved. The signal-to-noise ratio of the spectrum 
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Fig. 4. Contour maps of the observed CS (2-1), CS (3-2), 
H 13 CO+(l - 0) andDCO+(2 - 1) distributions in CG ^su- 
perposed on a grey scale map of the C 18 (2-1) emis- 
sion. The positions observed in each molecule are in- 
dicated. The corresponding SEST beam size at the ob- 
served frequencies is indicated in the upper left corner 
of each panel. In the CS (2-1) and CS (3-2) maps the 
lowest contour value is 0.5 Kkms -1 and the increments 
are 0.5 Kkms -1 and 0.3 Kkrns -1 , respectively. In the 
H 13 CO+(l - 0) and DCO+(2 - 1) maps the lowest contour 
value and its increment are 0.15 Kkms -1 . The maximum 
values of the CS (2-1), CS (3-2), H 13 CO+(l - 0) and 
DCO+ (2 - 1) emission are 2.2 Kkms -1 , 1.3 Kkms _1 , 
0.58 Kkms -1 and 0.77 Kkms -1 , respectively. The offset 
from the map zero position is shown on the axes. 



is not sufficient to make a meaningful fit of the line total optical 
depth, T tot . The estimated r to t upper limit in this position is 5. 
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Fig. 6. CS (2-1), (3-2) and C 34 S (1-0) spectra observed in the 
direction of CG 12-N (the two upper panels) and CG 12-S (the 
lower panel) 



3.3.2. CG 12-N 

In CG 12-N (position -40",160") the C ls O lines are nearly 
symmetric but not gaussian. The emission from the high 
density tracers is weak when compared to that observed in 
CG 12-S. 

The C 34 S line was observed in two positions in 
CG 12-N (Fig.®. The CS(2-1) and (3-2) line profiles in these 
positions are not gaussian and may be composed of two or three 
components. The peak intensity of the C 34 S (2-1) line is red- 
shifted with respect to the peak of the CS (3-2) line which itself 
is redshifted with respect to the CS (2-1) line. The latter line 
peaks approximately at the same velocity as the C 18 lines. 

The observed ionic lines (DCO+ (2-1), H 13 CO+ (1-0) and 
N2H + (1-0)) are offset +0.5 kms -1 in velocity with respect 
to the C 18 lines and peak at a velocity where especially the 
C ls O (2-1) line intensity is low. 

N2H + (1-0) was observed in the same two positions as 
C 34 S. Contrary to the CS lines, the observed N2H + lines are 
identical within the noise (Fig. [7] two upper panels). The SEST 
beam is the same at CS (2-1) and N2H + (1-0) frequencies. 
A hyperfine component fit to the CG 12-N N2H+ lines gives a 
0.65 kins -1 broad line at —5.41 kms" 1 . The estimated r to t 
upper limit in position (—40",— 160") is 2. 

To rule out that the observed intensity /velocity structure de- 
scribed above is due to calibration, pointing or frequency set- 
ting problems, the observations were carefully checked by ob- 
serving the C 18 and the DCO+, H 13 CO+ and CS lines after 
each other and making pointing checks before and after obser- 
vations. Also, e.g., the line pair DCO+ (2-1) and C 18 (1-0) 
could be observed simultaneously with a dual receiver, thus ex- 
cluding relative pointing errors. 
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Fig. 7. N2H" 1 " (1-0) spectra observed in the direction of 
CG 12-N (the two upper panels) and CG 12-S (the lower 
panel). The seven component hyperfine structure fit and the re- 
sulting residuals are shown 

4. Comparison with observations at other 
wavelengths 

The C 18 (2-1) contour map superposed on the blue SERC- 
J DSS image is shown in Fig. [8] The bright optical re- 
flection nebula NGC 5367 lies in front of CG 12-S whereas 
CG 12-N is in the direction of an optically heavily obscured 
region which is clearly seen in the inset. A 1.2 mm source 
(Haikala 2006, in preparation) coincides with the position of 
the CG 12-N C 18 maximum. 

A Ks band image (Haikala 2006, in preparation) of CG 12 
is shown in Fig. [9] Superposed on the image are the contours of 
the 12 CO (2-1) molecular outflow line wing emission (White 
1993). In the insets the C 18 (2-1) and DCO+ contour maps 
are superposed together with the positional uncertainty ellipses 
of IRAS point sources IRAS 13546-3941 and IRAS 13547- 
3944. IRAS 13547-3944 does not coincide with the binary 
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Fig. 5. Spectra of C 18 and molecules sensitive to number density towards the same positions as shown in Fig. [3] Only the 
unblended hyperfine N2H + (1-0) (Fi, F= 1,0 — > 1,2) component is shown. The velocity tick indicates a velocity of —6.5 
kms -1 . C 18 (1-0), CS (2-1) and H 13 CO + lines are plotted using continuous line. The panel lower right is a blowout of the 
upper right panel showing only the lines tracing high number density. 



h4636 (the base of the binary star emission is saturated due 
to the intensity scale which was chosen to emphasize the low 
intensity surface emission). 

A cone like nebulosity with a bright head is seen projected 
on the DCO + core in Fig. |9] The apex of the cone is located 
just off the tip of the red shifted and below the end of the blue 
shifted, collimated CO outflow lobes of White (1993). A com- 
pact 1.2 mm source (Haikala 2006, in preparation) coincides 
with the DCO + core. The core does not have an associated 
point source but this could be due to the low spatial resolution 
of the IRAS satellite. A faint source would be masked by the 
strong nearby source IRAS 13547-3944. The NIR cone could 
be associated with the driving source of the outflow. The centre 
of the molecular outflow is offset from the position of IRAS 
13547-3944 by ^20". It is unlikely that the point source is as- 
sociated with the large outflow but it could be the driving force 
of the strong redshifted outflow lobe located ~1' NW of the 
point source nominal position. 

Small areas of the size of the SEST HPBW at the C 18 (1- 
0) frequency in the centres of CG 12-N and CG 12-S have 
been mapped in the C 18 (3-2) (Haikala et al. 2006). The 
HPBW of the C 18 (3-2) observations is 19" which is simi- 



lar to the SEST HPBW of 24" at the C 18 (2-1) frequency. 
In CG 12-S the C 18 (3-2) emission has a strong maximum in 
the same position as the maxima in the lower C 18 transitions. 
The line profile in position (—20",— 20") agrees with the 
C 18 (2-1) line (Fig.0. The C 18 (3-2) T MB peak line tem- 
perature in this position is, however, 10 K which is nearly 
twice the C 18 (2-1) T MB of 5.7 K. In CG 12-N (position 
(-40", 160") the C 18 (3-2) line profile is similar to that of 
C 18 (1-0) but the Tmb peak line temperature is 3 K which 
is lower than the corresponding values observed in C 18 (1-0) 
and (2-1) which are 4.3 K and 3.8 K, respectively. 



5. C 18 fine scale structure 

The half widths of the observed C ls O lines in CG 12 cores 
range from 1.0 kms" 1 (CG 12-S) to 1.3 kms" 1 (CG 12-N). 
The C 18 channel-maps (Fig. 1111 and observations of other 
molecules indicate, that the lines consist of more than one ve- 
locity component. The division of the cloud according to its ap- 
pearance in the C 18 line integral maps (Fig.Q into only three 
components, CG 12-N, CG 12-S and CG 12-SW, is therefore 
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Fig. 8. Contour map of the C 18 (2-1) emission superposed 
on the optical (DSS blue) image of NGC 5367. The image is 
also shown in the inset but using an intensity scale which better 
brings out the surface emission. The linear size of the image is 
32'by 32'. 




1 0-1-2 
K.A. offset (arcmin) 



Fig. 9. Ks-band image of CG-12 (Haikala 2006, in prepa- 
ration). Contours of the blue (white contours) and red- 
shifted (black) 12 CO (2-1) line wing emission (White 1993) 
are also shown. The two insets show the CG 12-S and 
CG 12-N regions in detail. Superposed are the contours of 
C 18 (2-1) (black) and DCO+ (white) emission. The po- 
sitional uncertainty ellipses of the IRAS point sources IRAS 
1 3547-3944 in CG 12-S and IRAS 13546-3941 in CG 12-N are 
also shown. 



too coarse. The fine structure within the individual maxima is 
not seen because the lines are heavily blended in velocity. 

5. 1. Positive Matrix Factorization 

Positive Matrix Factorization (PMF) has been used by 
|Juvela et al. (1996)| and |Russeil at al. (2003)| in the analysis of 
molecular line spectral maps. PMF assumes, that the ensem- 
ble of input spectra are composed of n individual line com- 
ponents (factors). Unlike the Principal Component Analysis 
PMF assumes that the individual components are positive. This 
makes the interpretation of the results more straightforward 
than in Principal Component Analysis where the results may 
contain negative components. No other assumptions are made 
of the shape of the factors. Each input spectrum can be recon- 
structed from the n PMF output factors by multiplying each 
factor by the weight assigned to it by PMF at the position and 
adding up the multiplied factors. A comparison of the results 
obtained with PMF analysis with those obtained using the anal- 
ysis of channel maps and by the use of Principal component 
analysis is given |Russeil at al. (2003)| 

PMF has been applied separately to the north and south 
parts of CG 12 but there is a small positional overlap between 
the two areas used in the analysis. The same dataset, which is 
shown in Fig.^Jwas used in the analysis, i.e, the C 18 (2-1) 
data is binned to the same channel width as the C ls O (1-0) 
data. The PMF analysis is presented in AppendixlAl 



6. Discussion 

6. 1. C 18 as a tracer of molecular material 

C ls O emission is generally considered a good tracer of 
large scale structure of dark clouds and globules. However, 
C 18 photodissociation in the cloud envelopes and molecular 
depletion in the dense and cold cloud cores may restrict the 
number density interval where C 18 emission can be used as 
a direct measure of H2 column density (especially if the LTE 
approximation is used). Unlike 12 CO, the C 18 molecules are 
not shielded against photodissociation by strong H2 lines and 
further, the C 18 self-shielding is weaker than that of 12 CO. 
For C 18 the self-shielding is most efficient for the two low- 
est rotational levels with the largest populations. According to 
[Warin et al. (1996)1 the consequence is that the J=l-0 transi- 
tion becomes thermalized, whereas the higher transitions re- 
main subthermally excited in the cloud envelopes. In the other 
extreme (pre-stellar cloud cores and protostellar envelopes) 
the CO molecule may vanish from gas phase because of de- 
pletion onto dust grains. Further complications in interpreting 
C 18 data are the possible strong gradients in the CO excita- 
tion temperature on the line of sight due to heating of the gas 
by newly born stars and protostellar objects. Time and tem- 
perature dependent chemistry like, e.g., deuterium fractiona- 
tion, also complicates the comparison of C 18 data with that 
of other molecular species. 

Depletion of molecules on dust grains takes place in the 
dense, quiescent and cold molecular cloud cores. CO and CS 
are among the first molecules to disappear from the gas phase 
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but nitrogen bearing species like, e.g., N 2 H + and ammonia, 
seem to be able to resist depletion ( Tafa lla et al. 20 02 ). The 
N 2 H+and DCO+ lines in CG12-Npeak at a velocity of 
5.4 kms -1 where the C 18 line emission is weak (Sect. EOl . 
It is likely, that strong CO depletion has taken place in the gas 
traced by these two molecules in CG 12-N . 

In CG 12-S the C 18 maximum is spatially offset from the 
DCO + core. Deuterium fractionation reactions are favored in 
cold gas (e.g., Herbst 1982 1 and therefore the relative abun- 
dance of DCO + is enhanced in cold cloud cores. The relatively 
strong DCO + emission observed in CG 12-S, the DCO + core, 
could indicate that fractionation has indeed taken place. If the 
molecules in CG 12-S were heavily depleted in the cold DCO + 
core one would not expect CS (3-2) or H 13 CO + to outline the 
core like they do (also DCO + will be finally depleted). 

The observed 12 CO (1-0) and (2-1) line temperatures 
in CG 12-S (Fig. [5} are high. The lines are self reversed 
so the actual line peak temperatures must be even higher. 
This indicates that the CO excitation temperature is in ex- 
cess of 30 K in the part of the cloud which is traced by 
the observed 12 CO emission. White (1993) suggests that ei- 
ther h4636 or IRAS 13547-3944 is the heating source. The 
12 CO and 13 CO optical depth is high and unlike C 18 these 
CO isotopoloques are likely to trace only the surface of 
the molecular cloud associated with CG 12. The observed 
C 18 (3-2) T M b peak temperature of 11 K in CG 12-S also 
points at a high C 18 excitation temperature (Haikala et al. 
2006). However, the observed relative C 18 (2-1) and (1-0) 
line intensities in CG 12-S are not compatible with C 18 T ex 
values higher than 20 K. 

6.2. PMF: The interpretation 

One should be cautious in interpreting the PMF results. Even 
though the fit to the C 18 (1-0) and (2-1) data is good, all 
the PMF factors do not necessarily describe real cloud compo- 
nents. The structure of the cloud can, however, be discussed 
with some confidence when the PMF fit results are consid- 
ered with the information provided by other available molec- 
ular line, mm continuum and NIR-FIR data. 

6.2.1. CG 12-S 

The location of the CG 12-S C 18 (2-1) PMF factor maxima 
relative to the NIR cone (and the mm continuum source) in the 
centre of the DCO + core is shown in Fig. ^3 From East to 
West the factors are Is, 2s and 3s. The centres of the maxima 
differ in position by about one SEST beam size at 220 GHz, 
24". 

The data would seem to implicate that CG 12-S is frag- 
mented into three individual cores. If this were the case one 
would expect that the 2s factor would coincide with the DCO + 
core because the 2s centre of line velocity is the same as that of 
DCO + . Even though the 3s factor is the strongest in CG 12-S it 
is detected only as an asymmetry in the DCO + and CS lines in 
Figs.|5]and[6] CO is ubiquitous and easily excited and therefore 
it can be detected at much lower densities than the high density 
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Fig. 10. The CG 12-S C 18 (2-1) Is, 2s and 3s (from East to 
West) factor peak emission positions (0.8 K kms -1 contours) 
and DCO + core position (white) superposed on the SOFI Ks 
band image. Offsets from the map zero position are shown in 
arcminutes. The SEST beam at the C 18 (2-1) frequency is 
0.37 arcminutes 



tracers, DCO+ and CS (3-2), which trace the DCO+ core. The 
bulk of the observed C 18 emission may have its origin in the 
part of the cloud where the density is lower and the excitation 
temperature higher than in the DCO + core. 

This leads to a model where most of the observed 
C 18 emission traces rather the surface of the dense and cold 
cloud core (the DCO + core) or a separate, adjacent cloud com- 
ponent, than the core itself. The observed CS (2-1) distribution, 
which is similar to C 18 in CG 12-S (Fig.|4ji, is in accord with 
this model. Much of the observed CS (2-1) emission is known 
to originate in cloud envelopes and not in the dense cores. The 
CS (3-2) effective critical density is higher than that of CS (2- 
1) and therefore it traces the dense gas deeper in the cloud than 
CS (2-1). 

The possible interaction of the collimated molecular out- 
flow detected by |White (1993)| with the cloud core further com- 
plicates the interpretation of the molecular line data. The centre 
of the outflow is located in the direction of the DCO + core and 
it is not unlikely that the outflow can produce weak C 18 line 
wings and that it can raise the CO excitation temperature lo- 
cally. Weak line wings are also observed in the CS (2-1) line in 
Fig. |3 It is argued that the weak C 18 red and blue shifted line 
wing emission near the DCO + core is produced by the interac- 
tion of the collimated molecular outflow with the parent cloud. 
The line wings would be described by the PMF factors Is and 
4s. For the latter factor this would consists only of the emission 
immediately to the west of the DCO + core. 

The CG 12-S C 18 (3-2) map in Haikala et al. (2006) cov- 
ers only the very centre of CG 12-S. The C 18 (3-2) line peak 
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velocity and the intensity distribution is in accordance with the 
elongated shape of the 3s component in Fig. [TO] The maximum 
of the 3s factor is compact and lies on the outflow axis only 
r^2Q" away from the apex of the NIR cone (Fig.|9jl, the putative 
position of the outflow driving source. This suggest, that the 
high 3s C 18 line intensities and the molecular outflow could 
be connected. 

6.2.2. CG 12-N 

PMF produces a seemingly straightforward solution for 
CG 12-N. However, the available data from molecules other 
than C 18 and the 1.2 mm continuum emission deny such a 
simple solution. As argued in Sect. 16. H it is likely that the cloud 
component traced by DCO + and N 2 H + is heavily depleted. 
If this is the case C 18 only probes the undepleted part of 
CG 12-N. 

6. 3. Cloud physical properties 

A straightforward division of the CG 12 molecular cloud 
into three homogeneous components, CG 12-N, CG 12-S and 
CG 12-SW, is not possible. The good velocity resolution and 
the high signal to noise ratio of the data allows one to divide the 
C 18 data into separate components. However, the derivation 
of the cloud or core physical properties calls for a detailed three 
dimensional non-LTE model which takes into account both the 
molecular depletion, the varying density and excitation con- 
ditions. Such a model is beyond the scope of this paper and 
will be left for the future. The LTE-approximation approach is 
used instead of a sophisticated model to make a zeroth order 
estimate of the masses of the three C 18 maxima, CG 12-N, 
CG 12-S and CG 12-SW. An average C 18 excitation temper- 
ature for each maximum is estimated using the observed rela- 
tive C 18 (2-1) and (1-0) line intensities. 

6.3.1 . Mass and column density estimation 

The observed TJ line temperatures were converted into main 
beam temperatures using the beam efficiencies in Table 
^ The C 18 O excitation temperatures were estimated from 
the C 18 (1-0) and (2-1) data assuming optically thin 
C 18 emission and LTE. This assumes that the observed emis- 
sion in both transitions originates in the same volume of gas 
at a constant excitation temperature. The observed C 18 line 
T MB (2-Y)I T AI b(1-0) ratio is 1.8 in CG 12-S, compatible 
with an excitation temperature near 15 K. In CG 12-N and 
CG 12-SW this ratio is one or less indicating excitation tem- 
peratures of the order of 10 K. 

The masses of CG 12-N, CG 12-S and CG 12-SW are 
calculated without dividing them into smaller components. 
The spectra between declination offsets (including the lim- 
its) — 120" and +60" were assigned to CG 12-S. The spectra 
above and below these limits were assigned to CG 12-N and 
CG 12-SW, respectively. An average C 18 excitation temper- 
ature of 10 K is now assumed for CG 12-N and CG 12-SW and 
15 K for CG 12-S. The adopted C 18 abundance is 2.0 10~ 7 . 



This results in calculated total masses of 34, 96, and 110 
Mq for CG 12-SW, CG 12-S and CG 12-N, respectively, when 
C 18 (1-0) data is used. If C 18 (2-1) data is used the masses 
are 14, 49 and 53 Mq. The maximum Tmb line integrals can 
be used to estimate the H 2 column densities in the correspond- 
ing beams. The calculated maximum column densities in the 
C 18 (1-0) beam are 1.3 10 22 , 1.8 10 22 and 2.7 1 22 cm" 2 
for CG 12-SW, CG 12-S and CG 12-N, respectively. For the 
C 18 (2-1) beam the corresponding values are 6.4 10 21 , 1.6 
10 22 and 1.8 10 22 cm- 2 . 

The masses calculated from the C 18 (1-0) data are ap- 
proximately twice higher than from the (2-1) data. The LTE 
approximation assumes that all the C 18 rotational states are 
thermalized. However, according to Warin et al. (1996) only 
the lowest C ls O rotational state is thermalized in dark clouds 
and the higher states are subthermally excited. The subthermal 
excitation would be stronger in the relatively low density cloud 
envelope than in the more dense core. The subthermal level 
population leads to an underestimation of the cloud mass when 
calculated from the C 18 (2-1) data. This could be a partial 
explanation for the large discrepancy between the masses cal- 
culated using C 18 (1-0) and (2-1) data. 

Strong C 18 (3-2) emission (maximum Tmb ~H K) 
distributed similar to the 3s PMF factor was detected by 
IHaikala et al. (2006)| in the centre of CG 12-S . The C 18 (3- 
2) data was modelled with a compact (60" to 80" diam- 
eter) and hot (80 K <T ex < 200 K) optically thin clump 
of ~ 1.6 Mq. The high temperature was derived from the 
T MB (C 18 (3-2))/T MS (C 18 (2-1)) ratio. However, T MB 
(C 18 (2-l))/T MB (C 18 (1-0)) ratio for the CG 12-S PMF 
factor 3s is not compatible with such a high temperature. This 
could be due to, e.g., contribution from a subthermally ex- 
cited cool cloud envelope to the C 18 (1-0) emission. The LTE 
mass calculated from the C ls O (2-1) 3s component within the 
area modelled in Haik ala et al. (2006)1 is 4.1 M Q when a T ex 
of 15 K is assumed. The discrepancy in the calculated ex- 
citation temperatures and masses highlights the uncertainties 
when LTE approximation is used and only data of the two low- 
est C 18 transitions are available. In contrast to CG 12-S the 
C 18 T ex and the LTE mass of CG 12-N correspond well to 
the modelled values using the the three transitions (Haikala et 
al 2006). 



6.4. Summary 

The analysis of the spectral lines presented above deliver a 
complicated picture of CG 12. Even though the spectral lines 
are relatively narrow the analysis reveals a rich structure in line 
shape and velocity. Probable depletion of molecules on dust 
grains in CG 12-N further complicates the interpretation. The 
cloud parameters derived using the LTE analysis are highly un- 
certain. However, if the molecular line observations are com- 
bined with the available optical, NIR, FIR and mm continuum 
data the cloud structure can be discussed with some confidence. 
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6.4.1. CG 12-N 

CG 12-N harbors a compact, cold mm source and the detected 
relatively weak molecular emission in high density tracers is 
probably associated with this source. Much of the molecular 
material associated with this cloud component is likely to be 
highly depleted onto dust. The observed strong C 18 emission 
towards CG 12-N originates therefore in the envelope of this 
depleted core or in a separate entity seen in the same line 
of sight. 12 CO observed in the direction of CG 12-N has line 
wing emission indicating molecular outflow. It is not known if 
this outflow is connected to the outflow in CG 12-S or if it is 
local and originates in CG 12-N . 

6.4.2. CG 12-S 

The bulk of the observed C 18 emission does not trace the gas 
associated with the compact cloud core detected in H 13 CO + , 
DCO+ and CS (3-2) lines. Most likely the C 18 emission 
traces only the surface of this core. The moderate C 18 line 
wing emission, probably due to interaction of the highly col- 
limated molecular outflow with the surrounding gas, further 
complicates the interpretation. 

The molecular line data presented in this paper combined 
with the molecular outflow data from |White (1993)| the NIR 
imaging and mm dust continuum data (Haikala 2006, in prepa- 
ration) shows that the outflow centre coincides with the mm 
continuum source and the NIR cone in the centre of the DCO + 
core. The strong point source IRAS 13547-3944 is offset from 
the mm source and from h4636. This point source could, how- 
ever, be the driving source for the strong redshifted outflow 
lobe north of CG 12-S. 

6.4.3. CG 12-SW 

CG12-SW is inconspicuous when compared to the two stronger 
C 18 maxima. No signs of star formation have been found in 
its direction and it is therefore either still in pre-star forma- 
tion phase or its density-temperature structure is such that no 
star formation will take place. The arc like feature which con- 
nects CG 12-S and CG 12-SW at velocities from -6.6 kms" 1 
to -6.2 kms -1 in Fig. [TT] suggests that CG 12-SW might be 
connected with CG 12-S. 



7. Conclusions 

We have performed a detailed, high signal-to-noise ratio, mm 
line study of CG 12 in various molecular transitions, princi- 
pally of C ls O (1-0) and C ls O (2-1), as well as in molecular 
lines probing dense material, and have obtained the following 
results: 

1. The C ls O line emission is distributed in a 10' North- 
South elongated lane with two strong, compact maxima, 
CG 12-N and CG 12-S, and a weaker maximum, CG 12-SW. 

2. High density tracers CS (2-1), (3-2), DCO+ and 
H 13 CO + are detected in both strong C 18 maxima. Emission 
from these molecules is weak in CG 12-N but in CG 12-S it 



defines a compact core (referred to as DCO + core) which is 
spatially offset from the C 18 maximum. 

3. The emission from the high density tracers in 
CG 12-N takes place at a velocity where emission from 
C 18 is weak. The molecules associated with the cloud com- 
ponent detected in high density tracers are likely to be heavily 
depleted. 

4. Positive Matrix Factorization was applied to study 
the cloud C 18 fine scale structure. The observed strong 
C 18 emission in CG 12-N (PMF factor 2n) originates in the 
envelope of the depleted cloud component or in a separate en- 
tity seen in the same line of sight. In CG 12-S the most intense 
PMF factor 3s traces warm gas on the surface of the DCO + 
core or a separate adjacent cloud component. 

5. The driving source of the collimated molecular outflow 
detected by |White (1993)| lies in the DCO+ core. 

6. The average C 18 OLTE mass is -80 M for 
CG 12-N, -70 M for CG 12-S and -20 M Q for CG 12-SW. 
These numbers can only be considered as a zeroth or- 
der estimate because of the uncertainty in defining the 
C 18 excitation temperature and possible molecular depletion 
in the C 18 maxima. 

7. If the distance to CG 12, 630 pc, is correct the linear 
size and the mass of this cometary globule approaches that of a 
typical low mass star forming region like eg. Chamaeleon I. 
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Appendix A: Positive Matrix Factorization 
analysis of CG 12 

A.1.CG 12-N 

The results of the PMF, assuming three factors, to observed 
C 18 (1-0) and (2-1) data in CG 12-N are shown in Figs.rO 
and lA.2l The uppermost panels show the three fitted PMF fac- 
tors. The area under each factor is normalized to one. The his- 
togram superposed on the rightmost factor in Fig lA.2l indicates 
the velocity resolution of the input data. The intensity distri- 
bution of the PMF factors in CG 12-N are shown in the lower 
three panels. 

The PMF reproduces well the velocity structure displayed 
in Fig. El The fits to the C 18 (1-0) and (2-1) data are inde- 
pendent of each other. Therefore it is encouraging that the PMF 
output factors for both C 18 transitions are similar in velocity 
and in spatial distribution. Because of the unambiguity of the 
fits for both transitions the factors, centered at —6.6, —5.9 and 
—5.4 kms -1 , are referred to in the following as In, 2n and 3n, 
respectively. 

The PMF factors, In, 2n and 3n, are distributed in a north- 
south oriented ridge, in a pear shaped body and a narrow 
east-west ridge, respectively. The 3n factor lies at the ve- 
locity where the emission from the H 13 CO + , DCO + , and 
N2H + is at maximum. Even though this factor is not read- 
ily evident in the individual spectra PMF finds it in both 
C 18 transitions. Factors 2n and 3n are symmetric but fac- 
tor In has blue and redshifted wing like emission, especially 
in the (1-0) transition. According to the PMF fit the high 
T% (C ls O (l-0))/TJ (C ls O (2-1)) ratio observed in the direc- 
tion of CG 12-N is due to the factor 2n. 

If the number of factors in the PMF is chosen to be four 
PMF divides in essence the factor In into two, leaving the two 
other factors untouched. The spatial distribution of the split up 
factors is similar to that of the In in the three-factor PMF. If the 
number of factors is further increased to five the result is similar 
to that with four factors plus a fifth factor which has nearly 
zero intensity, 'empty' field, in the map. The three factors are 
therefore sufficient to produce the velocity structure evident in 
Fig.^Jand increasing the number of factors does not improve 
the fit significantly. 

A.2. CG 12-SandCG 12-SW 

The results of the PMF fits, assuming four factors, to both ob- 
served C 18 transitions in CG 12-S are shown in Figs. lA.3l and 
IA.4I The location of the DCO + core is shown with a dashed 
contour in the figures. The factors, centered at —6.8, —6.4, 
—6.2 and —5.8 kms -1 , will be referred to as Is to As, respec- 
tively. 

The most redshifted factor 4s covers the very northern 
part of the region and is close in velocity to the 2n factor in 
CG 12-N. It is natural to consider that this factor is due to emis- 
sion extending from CG 12-N to CG 12-S. Also the very faint 
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Fig. A.l. Upper panel: the basic spectral profiles (factors) cal- 
culated by the PMF in CG 12-N (three factor fit), C 18 (1- 
0) data. Lower panels: maps of the weights of the basic fac- 
tors. The lowest contour level and increment are 0.2 K km s _1 . 
Offsets from the centre position of the original C 18 map are 
shown in arc minutes on the x and y axes. 



emission observed in the northern part of the field in other three 
factors is at least partly due to emission from CG 12-N. 

The emission from the most blueshifted factor Is is concen- 
trated just below the (0,0) position. The 2s factor represents the 
arc seen in Fig. which connects CG 12-S and CG 12-SW . 
There is a small size local maximum in the C 18 (2-1) 2s fac- 
tor north of the DCO + core. Even though this local maximum 
is not seen in the C 18 (1-0) 2s factor, there is extended emis- 
sion at its location. The southern part of the arc is more intense 
in the (1-0) transition than in the (2-1). The maximum of the 
southern extension coincides with CG 12-SW. The emission 
from the 3s factor peaks west of the DCO + core and is seen in 
both transitions. 
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Fig. 11. T'\ channel map of the C ls O (1-0) emission (upper three rows) and C ls O (2-1) (lower three rows). The colour scale is 
the same for both transitions. The LSR velocity is indicated in the upper left corners of the panels. Each pixel corresponds to a 
single observed position but the C ls O (2-1) data is binned in this figure to the same channel width as the (1-0) transition, 0. 1 16 
kms -1 . The highest intensity in the panels is ~ 3.0 K. The cross in the panels is located in the map centre position. 
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Fig. A.5. The PMF fit to the C 18 lines in the selected positions positions in Fig.|5] The observed C 18 (2-1) and (1-0) lines 
are plotted with a heavy histogram and line, respectively. The green and red lines show the fitted PMF C 18 (2-1) and (1-0) 
factors, respectively. The residuals of the fits are shown in the lower part of each panel. 



A. 3. Individual line profiles 

The PMF fits to the C 18 lines in the same selected positions 
as in Fig. [5] are shown in Fig. IA.5I The residuals obtained by 
subtracting the added up PMF factors from the observed lines 
are shown below the lines in the figure. A grey scale map of 
all the PMF fit residuals is shown in Fig. IA.6I One should 
stress that, unlike a gaussian multicomponent fit to a single 
spectrum, PMF fits, in this case 3 (CG 12-N) or 4 (CG 12-S) 
factors into all the input spectra simultaneously. Neither the 
shape nor the velocity of the factors is fixed in the fit. The two 
C 18 transitions were fit separately and PMF could have used 
factors differing in profile and velocity for the two transitions. 
The fitted PMF factors are, however, similar both in shape and 
velocity. The residuals shown in Figs. IA.51 and IA.61 are small 
and demonstrate that PMF produces a good fit to the data. 



The largest residual in Fig. IA.5I takes place in the 
blue shifted side of the C 18 (1-0) line in the centre of 
CG 12-N (position -40",160"). This would imply that the 
blue shifted wing of factor C 18 (1-0) In is too strong. A 
closer comparison of the input data and PMF results reveals 
that similar deviations take place for three other C 18 (1— 
0) spectra around the position above. The C 18 (1-0) line 



profiles farther away from the CG 12-N centre actually do 
have a more prominent blue shifted wing than the four above 
mentioned spectra. The large number of spectra with good 
fits outweighs the less optimal fit for the spectra in the cen- 
tre of CG 12-N. The residuals shown in Fig. IA.5l for the re- 
maining four spectra in CG 12-S are smaller than in posi- 
tion (-40", 160"). 
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Fig.A.2. As Fig.|AT]but for C 18 (2-1) data. Factor 3n in 
the upper panel has been plotted also as a histogram to indicate 
velocity resolution of the input/output spectra. 
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Fig.A.3. Upper panel: the basic spectral profiles calculated 
by the PMF in CG 12-S (four factor fit), C 18 (1-0) data. 
Lower panels: the maps of the intensity of the basic factors. The 
factor numbers correspond to the numbers in the upper panel. 
The lowest contour level and increment is 0.2 K km 
dashed contour outlines the DCO + core. 
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Fig. A.4. As Fig.lAJlbut for C 18 (2-1) data. 
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Fig.A.6. The PMF fit residuals in the K km s 1 scale. 



